Many-body effects on the Rashba-type spin splitting in bulk bismuth tellurohalides 
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We report on many-body corrections to one-electron energy spectra of bulk bismuth 
tellurohalides — materials that exhibit a giant Rashba-type spin splitting of the band-gap edge states. 
We show that the corrections obtained in the one-shot GW approximation noticeably modify the 
spin-orbit-induced spin splitting evaluated within density functional theory. We demonstrate that 
taking into account many-body effects is crucial to interpret the available experimental data. 

PACS numbers: 71.15.m, 71.20.b, 71.70.Ej 



I. INTRODUCTION 

Recently, bismuth tellurohalides have been discovered 
as a new very promising class of materials for semi- 
conductor spintronics applications li"— Due to inversion 
asymmetry of bulk crystal potential, these semiconduc- 
tors show a giant Rashba-type spin splitting of the bulk 
states. A study of bismuth tellurohalides has been trig- 
gered by Ref. [ll where the band structure of BiTel 
has been investigated with angle-resolved photoemission 
spectroscopy (ARPES). The spin-split state experimen- 
tally observed in Ref. [l| was attributed to the bulk con- 
duction band minimum (CBM) shifted downward due 
to band bending in the near-surface region and confined 
in the respective accumulation layer. The splitting of 
such a CBM was characterized by the momentum offset 
kn ^ 0.052 and the Rashba energy Er ~ 0.1 eV. 
These experimental values for the Rashba parameters 
were found to be in close agreement with those obtained 
from ah initio density functional theory (DFT) calcula- 
tions of the bulk BiTel band structure (0.05 and 
0.113 eV, respectively, as evaluated in Ref. [ll). This fact 
was interpreted in favor of a bulk nature of the observed 
state. 

In a DFT study i of bulk BiTeCl, BiTeBr, and BiTel, 
it was revealed that these polar layered semiconductors 
all exhibit similarly large Rashba-type spin splittings. 
Additionally, in Ref. and later in Refs. i and H it 
was shown that a free-electron-like spin-orbit-split sur- 
face state emerges in the bulk band gap at the (0001) sur- 
face of all mentioned semiconductors (hereafter we mean 
the Te-terminated surface). Similarly to the bulk CBM, 
this surface state possesses a Rashba-type spin split- 
ting and the corresponding spin texture. Quantitatively, 
e.g., in the case of BiTel, the surface-state spin split- 
ting is described theoretically by the Rashba parameters^ 
Er = 0.121 eV and = 0.068 1"^ The existence of 



this surface state has been experimentally corroborated 
in Refs.[^and0. Additionally, in these works it was ques- 
tioned if the state observed by the ARPES in Ref. [l| was 
really a bulk-derived state confined in the band-bending 
accumulation layer (see also Ref. @). However, when 
comparing the original ARPES measurements with the 
above theoretical Rashba parameters for the surface state 
we find that they disagree. As to BiTeBr and BiTeCl, 
their surface electronic structure has been investigated by 
ARPES very recently in Ref. H. Analogously to BiTel, 
signatures of spin-split states of a Rashba-type nature 
have been observed and interpreted in a similar way as 
in Ref. [l]. For these states, unlike BiTel, experimental 
values of the Rashba parameters reported in Ref. H can- 
not be supported well by the available theoretical data 
obtained for the bulk conduction-band bottom. 

When comparing the experimental ARPES data to 
band structures calculated within a Kohn-Sham (KS) 
DFT scheme one should keep in mind that, strictly speak- 
ing, the KS energies cannot be interpreted as excita- 
tion energies. For example, the fact that the KS band 
gap underestimates the experimental one often by a fac- 
tor of two is known as the hand-gap problem of DFT. 
Also, the dispersion of its edges can be affected. An ef- 
fective single-particle (quasiparticle) description of the 
many-body system can be retained if one employs an 
electronic self-energy that embodies the many-body ex- 
change and correlation effects. Already on the level of 
the GW approximation for the self-energy, calculations 
of many-body corrections to KS bands allow one to sig- 
nificantly improve theoretical results (see, e.g., Ref. [lO| ). 
However, ah initio GW calculations aimed at a study of 
the surface electronic structure of materials, where spin- 
orbit interaction plays an important role, are not feasible 
so far. 

In this paper, we present a theoretical ah initio study 
of the bulk band structure of BiTeX (X=I, Br, CI) tak- 
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ing into account many-body corrections evaluated within 
the GW approximation. We also distinguish between the 
local-density (LDA) and the generalized gradient approx- 
imation (GGA) for the exchange-correlation functional 
used in the calculation of the KS reference system. First, 
the quasiparticle corrections to the KS bulk states will 
provide an answer to the question of whether the cor- 
rected values of the bulk CBM Rashba parameters still 
match the values coming from the experimental observa- 
tion. Second, the fact that the surface-state properties 
are predetermined by the properties of the bulk band-gap 
edges^ allows one to estimate the effect of many-body 
corrections on the surface-state dispersion on the basis 
of the GW results for the bulk quasiparticle spectra. 

As a first result, we find that the GW approximation 
provides a more consistent picture than KS-DFT: While 
the Rashba parameters differ substantially for LDA and 
GGA, they are brought back into good agreement by 
the quasiparticle correction, demonstrating a weak de- 
pendence on the reference KS system that is used as the 
starting point for the GW calculation. 

We show that the quasiparticle corrections have a sig- 
nificant effect on the band gap and the Rashba parame- 
ters which determine the spin splitting of the band-gap 
edge states. For BiTel, the resulting GW values of the 
bulk parameters are no longer in close agreement to the 
aforementioned ARPES measurements of Ref. [l|. We sug- 
gest an estimate in order to predict changes in the dis- 
persion of the surface states induced by the many-body 
corrections. This estimate brings the theoretical and ex- 
perimental data back into good agreement. In the case 
of BiTeCl, the GW Rashba parameters are much closer 
to the experimental values than the KS results. There 
is only little difference between the GW bulk parame- 
ters and the estimate for the surface state, and both fall 
in the error range of the experiment. The quasiparticle 
correction to BiTeBr yields a mixed success. While it 
improves on the bulk parameter an, reflecting the ra- 
tio of the Rashba energy Er and the momentum offset 
knr with respect to LDA and GGA, the parameters Er 
and kji themselves are both underestimated, and the esti- 
mate for the surface state is no improvement in this case. 
Yet, the theoretical LBA+GW and GGA+GW values 
are again in very good agreement, while the correspond- 
ing LDA and GGA show notable differences. 



II. COMPUTATIONAL DETAILS 

We employ the full-potential linearized augmented 
plane wave (FLAPW) method as implemented in the 
FLEUR codeii within both the LDA of Ref. [H and 
the GGA of Ref. [l^ for the exchange-correlation (XC) 
functional. We use these two approximations in or- 
der to gain insight into the effect of different reference 
one-particle band structures on results obtained within 
the one-shot GW approximation realized by the SPEX 
code-ii The DFT calculations were carried out with the 



use of a plane- wave cutoff of k^ax = 4.0 bohr"-'^, an an- 
gular momentum cutoff of Imax = 12, and a 7 x 7 x 7 
F-centered fc-point sampling of the Brillouin zone (BZ). 
In order to treat quite shallow semi-core c?-states, the 
FLAPW basis was extended by inclusion of conventional 
local orbitalsji^ To more accurately describe high-lying 
unoccupied stateSfi^^ one local orbital per angular mo- 
mentum up to / = 3 was included for each atom. The 
linearization energy parameters for these local orbitals 
were set at the "center of gravity" of the respective par- 
tial density of states for each angular momentum I within 
the energy interval from the Fermi level up to a higher 
energy of 90 eV. In order to document the quality of con- 
vergence, we report GW results obtained with different 
numbers of states (Ni,) involved in the calculation of the 
dielectric matrix and the Green function and with two 
F-centered Monkhorst-Pack grids ( A''^ x Nk x Nk , where 
Nk = 5 or 6). The dielectric matrix evaluated within 
the random-phase approximation (RPA) was represented 
with the use of the mixed product basiSfi^ where we chose 
an angular momentum cutoff in the muffin-tin spheres of 
4 and a linear momentum cutoff of 3.5 bohr^^. The spin- 
orbit interaction was included into the GW calculations 
already at the level of the reference system. Not only the 
Green function but also the self-energy then acquire spin- 
off diagonal terms that lead to a many-body renormal- 
ization of the spin-orbit coupling, a phenomenon which 
is absent otherwise. 

We consider all bismuth tellurohalides in the hexagonal 
crystal structure with the lattice parameters reported in 
Ref. 1 1 3 and the atomic positions obtained theoretically in 
Ref. |3| by structural optimization. The crystal structure 
of the compounds is built up of alternating hexagonal 
layers Te-Bi-X stacked along the c axis. Each three layers 
form a three-layer (TL) block, and the distance between 
these blocks is about one and a half times greater than 
interlayer distances within the block. BiTeBr is supposed 
to be in an ordered phase with the atomic order similar 
to that in BiTel but with Br instead of I (see also Ref. 0)- 
The structure of BiTeCl differs from that of iodide in the 
sense that the stacking order of CI and Te layers alternate 
along the c direction resulting in a doubling of the unit 
ceU^ 



III. RESULTS AND DISCUSSION 

The obtained LDA bulk band structures are presented 
in Fig. [T] (upper row). As seen in the figure, within the 
H ~ A - L plane (or A" - F - M in BiTeCl, where the 
Brillouin-zonc is folded along F — ^ owing to the dou- 
bling of the unit cell) that is normal to the F — ^ line 
the bands are generally spin-split. Along the T ~ A line, 
bands are spin degenerate and show a relatively weak 
dispersion compared to that along the mentioned plane. 
This reflects a quasi-two-dimensional character of the 
corresponding states, which is caused by the pronounced 
TL crystal structure of BiTeX. In the case of BiTel and 
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FIG. 1: (Color online) Upper row: Bulk band structure of bismuth tellurohalides, as obtained within Kohn-Sham DFT with 
the use of different approximations (LDA and GGA) for the exchange-correlation functional. Shaded areas indicate the energy- 
momentum regions presented in the lower row. Lower row: The respective quasiparticle spectra evaluated within the GW 
approximation with the convergence parameters Nk = 5 and A^i, = 350 (for BiTeCl the number of bands is Nb — 580). 



BiTcBr, the band gap is formed by the minima of the 
conduction band and the maxima of the valence band in 
the vicinity of the A point of the Brillouin zone (BZ). 
The parameters that characterize the Rashba-type spin 
splitting of the conduction band in the A — L direction of 
the BZ in BiTel and BiTeBr are listed in Tables |T] and HIl 
respectively. These parameters, and also the band gap 
Egap, vary notably with the approximation chosen for the 
XC functional. The LDA gives systematically larger kn, 
Er, and an = 2Eji/kji and a smaller Egap than those 
obtained with the use of the GGA (note that the inverse 
relation between and Egap was revealed in Ref. 0). 
For BiTel, the LDA and GGA values of ka are quite 
close to the experimental values found in Refs. [H and H. 
Good agreement with experiment can also be found for 
the Rashba energy En evaluated within the GGA. As to 



the parameter an, which is quite sensitive to even small 
changes in kn and En, the GGA XC functional leads to a 
value that is in the error range of the experimental value 
4.3±0.9 determined from the measured momentum offset 
and Rashba energy in Ref. U For BiTeBr, the DFT cal- 
culation gives Rashba parameters of which only fc/j is in 
the error range of the corresponding experimental value 
(see TableHJ). In the case of BiTeCl (see TableHTII, where 
the extrema of the valence and conduction bands appear 
at the F point (see Fig. [T]), the Rashba parameters are 
again larger within the LDA than those obtained within 
the GGA. The latter are somewhat closer to experiment, 
while, in fact, only the values k^, and an fall within the 
error range of the experimental data. 

We performed calculations of the many-body correc- 
tions with different convergence parameters Nk and Nh- 
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TABLE I: The Rashba parameters (fc_R, Er, or) for the spin- 
orbit-spht conduction band and the band-gap width (-Egap) 
for BiTel, which were obtained without (LDA and GGA) 
and with (LDA+GW and GGA+GVK) many-body correc- 
tions evaluated within the GW approximation at different 
values of the convergence parameters A^^ and Ni,. An esti- 
mate of the Rashba parameters for the surface state (see the 
text) is also presented. The experimental values taken from 
Ref. are shown with the respective error range represented 
by numbers in parentheses. 



TABLE III: Same as in Table H] but for BiTeCl. 
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TABLE IL Same as in Table HI but for BiTeBr. 



Nk/Nb 


-E'gap, 


Er, 


kR, 


CtR, 




meV 


meV 


10-' A-i 


A-eV 






GGA 








310 


55 

GGA+GW 


34 


3.3 


5/350 


765 


20 


21 


1.9 


5/500 


760 


20 


21 


1.9 


6/350 


752 


20 

LDA 


21 


1.9 




207 


73 

LDA+GW 


37 


4.0 


5/350 


650 


22 


19 


2.2 


Estimate 




17 

Experiment 


27 


1.3 






42(10)" 


43(10)" 


2.0(7) 



"From Ref. H 



In Table IIVI we present the corresponding RPA dielec- 
tric constants eoo, which gives an indication of how the 
screening properties depend on the convergence parame- 
ters. As is clearly seen in the table, with the use of the 
GGA band structure an increase of the number of bands 
at a fixed A^^ has practically no effect on eoo- For BiTeBr, 
the dielectric constant does not show significant changes, 
when A^fc is increased at fixed number of bands. However, 
increasing at fixed Nf, has a stronger effect in BiTcI 
and BiTeCl: The dielectric constant grows by ^ 9% for 
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the former, while it decreases by ~ 5% for the latter. A 
similar situation is observed in the LDA case. An analysis 
has shown that for BiTel the LDA-based RPA calcula- 
tions of eoo with Nk = 7 and Nj, = 350 further increases 
the dielectric constant by ~ 5%. As we demonstrate be- 
low, the quasiparticle spectrum is much less sensitive to 
the changes in Nk and N^, and the fact that the dielec- 
tric constant is not fully converged with respect to the 
number of k-points (except for the BiTeBr case) does not 
affect the conclusions made from the obtained results. 

Now we discuss our results of the GW calculations, 
which are based on the LDA or GGA electronic states 
as a reference one-particle band structure. Fig. [T] (lower 
row) represents the GW results obtained with N^. — 5 
and Nb ^ 350 (580 in the case of BiTcCl). It is worth 
noting that the changes in Nk and Nh have hardly an 
effect on the quasiparticle spectra, which is reflected in 
Tables Hl lIIII As is clearly seen in the figure, the many- 
body corrections lead to a larger band gap, a smaller mo- 
mentum offset, and a steeper dispersion along the T — A 
line. Having a closer look at the values in the Tables 
Hl lIIII the first observation is that the many-body quasi- 
particle corrections increase the band gap considerably 
with respect to the LDA and GGA values, bringing it 
very close to the experimental value indeed. In the case 
of BiTel, the quasiparticle correction somewhat overesti- 
mates the experimental gap when the GGA reference is 
used, while the LDA+GW^ value is nearly on top of the 
experimental one. Second, the GGA+GW^ results listed 
in the tables demonstrate that the Rashba parameters 
arc well converged. That is why we have performed the 
LDA+GW calculations with A^^ = 5 and A^,, = 350 (580) 
only. Third, in the GW calculations the difference be- 
tween the Rashba parameters {Ef( and kn) obtained with 
the use of different approximations to the XC functional 
is smaller than in the DFT calculations. Thus, we can 
infer with confidence that taking many-body corrections 
to the DFT band structure into account leads to a re- 
duction of the Rashba energy and the momentum offset. 
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TABLE IV: The dielectric constant eoa {E _L c) obtained 
within the RPA for the considered semiconductors (for BiTeCl 
Nb is shown in parentheses). The experimental value of too 
for BiTel is taken from Ref. [IQ. 
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which characterizes the Rashba-type spin-splitting of the 
bulk band-gap edges. The ratio of these parameters, as 
represented by a^, also demonstrates a reduction except 
for BiTcI, where it becomes larger. 

In BiTel, the obtained GW results "worsen" the agree- 
ment with the experimental data. As a consequence, the 
interpretation done in Ref. [l| concerning the bulk nature 
of the Rashba-type spin-split state appears to be less le- 
gitimate in the light of the present GW values. On the 
other hand, in BiTeCl the GW corrections bring the theo- 
retical and experimental values of the Rashba parameters 
into better agreement. The case of BiTeBr is a peculiar 
one. Here, the GW values for the Rashba energy En 
and momentum offset fc^ underestimate the experimen- 
tal values by about the same factor so that their ratio, 
given by the parameter, is in very close agreement 
with the experimental data. 

Until the properties of quasiparticles in the surface 
state are calculated directly (which is a separate and very 
complicated problem at present), the presented results 
may serve as a basis for predicting the surface quasipar- 
ticle spectrum. Actually, as was experimentally estab- 
lished in Ref. the offset kn in BiTel found for the sur- 
face state is about 20% larger than that for the bulk con- 
duction band. The same change has been theoretically 
obtained in Ref. |3. Additionally, from DFT studies^i^ 
we know that the ratio between an for the surface state 
and that for the respective bulk state amounts to 0.66. 
(Note that for BiTeCl one gets practically the same value, 
whereas for BiTeBr it is equal to 0.56^) For consistency, 
a similar ratio for the Rashba energy En should be equal 
to 0.79. With these findings at hand, we can estimate 
Rashba parameters for the surface state at the BiTel 
(0001) surface. Assuming the LDA+GW values as a ba- 
sis, we enlarge kn by 20% and multiply En and an by 
0.79 and 0.66, respectively. As a result, we arrive at val- 
ues listed in Table U in the row marked as "Estimate" . 
This estimate restores the good agreement, which was 
lost upon taking the GW quasiparticle corrections into 
account. This can be interpreted in favor of a surface na- 
ture of the state experimentally observed in Refs. [l| and 
i 

Along with the dielectric constant and the band-gap 
width, the Rashba parameters estimated for the surface 
state in BiTel indicate that in spite of the fact that at 



the DFT level, as a rule, the GGA gives more accurate 
results for the band structure than the LDA does, in 
pairs with the one-shot GW approximation the latter 
becomes more preferable. Thus, with the use of the re- 
sults reported in Refs. [3, 0, and0 we can make a similar 
estimate for the surface states at the (0001) surface of 
BiTeBr and BiTeCl. As seen in Tables M and |llll the 
resulting kn and En values do not cause a change of the 
situation: in BiTeCl the modifications occur within the 
error range of the experimental values, and in BiTeBr the 
discrepancy with the experiment remains. This discrep- 
ancy could be caused by the fact that we use the lattice 
parameters a and c reported in Ref. [l^ for the disordered 
BiTeBr, where tellurium and bromine atoms randomly 
distributed within two layers facing Bi-atomic layer. An 
ordered BiTeBr has been grown very recently in Ref. |^ by 
the chemical vapor transport method. However, in Ref. |9| 
there is no information about the lattice parameters and 
atomic positions. We suppose that apart from surface 
effects a possible difference in the crystal structure can 
lead to the disagreement we observe. 



IV. CONCLUSION 

In conclusion, we have presented quasiparticle spec- 
tra of bulk bismuth tellurohalides. We have shown that 
taking into account many-body corrections to the KS 
band structure from DFT calculations leads to a no- 
table modification of the Rashba parameters {En, kn, 
and an), which characterize the spin-orbit splitting of 
the conduction-band minimum. We have revealed that 
the resulting GW values of these parameters are well con- 
verged at a moderate k-point sampling of the BZ and 
number of bands and demonstrate a weak dependence 
on the approximation chosen for the exchange-correlation 
functional. We have found that the GW corrections im- 
prove the bulk band-gap value considerably with respect 
to the KS value. It should be noted that this value turns 
out to be quite sensitive to the choice of the exchange- 
correlation functional on the level of KS DFT alone. We 
have shown that for BiTel the GW values of the Rashba 
parameters for the bulk conduction band worsen the good 
agreement with experiment that was reached at the DFT 
level. Thus, in the light of the GW results, the interpre- 
tation of the state observed by the ARPES in Ref. |l| to 
be of a bulk nature appears less likely. The suggested 
estimate of these parameters for the surface state has re- 
stored the good agreement, which may indicate a surface 
nature of the state in question. We have also found that 
in the case of BiTeCl the experimental data cannot be 
reproduced by DFT. However, both the GW values of 
the bulk Rashba parameters and estimates for those of 
the surface state match the experimental values within 
the error bars of the experiment. As to BiTeBr, neither 
the DFT nor the GW approximation and the further es- 
timate is able to reproduce simultaneously the available 
experimental data on En, kn, and an- 
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